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f-\ Abstract. We consider the entanglement between quantum field degrees of freedom inside 

and outside the horizon as a plausible source of black-hole entropy. We examine possible 
deviations of black hole entropy from area proportionality. We show that while the area law 
^^ holds when the field is in its ground state, a correction term proportional to a fractional power 
OQ of area results when the field is in a superposition of ground and excited states. We compare 
our results with the other approaches in the literature. 
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I— ' 1. Introduction 

I It is now well-known that the black-hole (commonly referred to as Bekenstein-Hawking) entropy 

^ is finite and is given by the relation [1] 
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^Bn = ^-^ ^here ^p = ^, (1) 



I A^ is the area of the black-hole horizon. It is infinite, if both the matter and gravity are 

t~^ purely classical! In other words, the finiteness of the black-hole entropy (and the validity of 

^^ second-law of thermodynamics) requires that the matter and(or) gravity to be quantized. 

J> Although the above relation has been derived by various approaches [2, 3, 4, 5], we do 

not yet understand: Why S'gjj is not extensive? What is the microscopic origin of black-hole 

entropy? Are there corrections to the Bekenstein-Hawking entropy? and How generic are these 

^ corrections? 

The purpose of this paper is an attempt to understand the generic corrections to the 
Bekenstein-Hawking entropy by assuming entanglement as a source of black-hole entropy. We 
show that while the area law holds when the field is in its ground state, a correction term 
proportional to a fractional power of area results when the field is in a superposition of ground 
and excited states. 
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The paper is organized as follows: In the next section, we briefly discuss various properties of 
entanglement entropy and provide motivation for the relevance of the entanglement entropy to 
S'gjj. In Sec. (3), we discuss the setup and the key steps in evaluating the entanglement entropy 
of scalar field propagating in flat space-time. In Sec. (4), we obtain the entanglement entropy 
for a superposed state and show explicitly that the power-law corrections become important for 
small horizon limit. We conclude with a summary in Sec. (5). In appendix, we show explicitly 
that at a fixed Lemaitre time, the Hamiltonian of a scalar field in a black hole background 
reduces to that in flat space-time. 

2. Entanglement entropy 

Let us consider a bipartite quantum system X = {12} and assume that X is in a pure state |^). 
The entanglement entropy {S^^J is defined as 

S^^,{^) ^ S{pi) = S{p2) (2) 

where 

S{pi) = -Tr[piln{pi)] = -'^Pnlnpn (3) 

n 

is the Von Neumann entropy, pn are the Eigen-values of pi obtained by solving the integral 
equation 

/•oo 

dx'pi{x, x')fn{x') = Pnfn{x) (4) 



and 

p^=Tr2mm) (5) 

is the reduced density matrix obtained by tracing over 2. Note that the entanglement entropy 
(2) is independent of the which degrees of freedom are traced over. 

In order to set the ideas, we consider a simple example — 2-coupled harmonic oscillator. 
The Hamiltonian of a coupled Harmonic oscillator is given by 
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P1+P2 + ko (xi + X2) + ki (xi - X2) (6) 



where {xi,pi) are coordinate and momenta of the oscillator, respectively, and ki is the 
interaction term which is assumed to be positive. The quantization of the above Hamiltonian 
is straight forward in the normal coordinates: 

(7) 
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uj^ = y%; uj^ = ^JkQ + 2ki; R^ = — < 1 . (9) 



Let us consider the situation in which we choose not to have any information about the 
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Figure 1. Plot of entanglement entropy of a coupled harmonic oscillator and R. 



oscillator 1. Mathematically, this corresponds to tracing over the oscillator 1, i. e. 

POO 

P2{X2,X2) = / dXi^Q{xi,X2)tpl{xi,x'2) 
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where 
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4(l + i?2) ' ' " 4(l + i?2) ' - VI + ^. 

Substituting (10) in the integral equation (4), the eigenfunctions and eigenvalues are given by 

„2^ 



fn{x) = Hn {^/ax) exp 



ax 



Pn = {l-OC 



Substituting the eigenvalues in (3), the entanglement entropy is given by 



5.„,(0 = -in(i-e) 
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ln(0 . 



(12) 
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In Fig. (2), we have plotted S^^^^. of the coupled oscillator in terms of R. It is interesting to 
note that the entanglement entropy vanishes when the oscillators are uncoupled {R= 1) while 
they diverge in the strongly coupled limit (i? = 0). 

But what has this digression to do with black-holes and what is the relevance of S^^^ to 
black-hole entropy? This can be understood from the fact that both entropies are associated 
with the existence of the horizons. Consider a scalar field on a background of a collapsing star. 
At early times, there is no horizon, and both the entropies are zero. However, once the horizon 
forms, 5gjj is non-zero, and if the scalar field degrees of freedom inside the horizon are traced 
over, S^^^ obtained from the reduced density matrix is non-zero as well. 



In the next section, we give the procedure for evaluating the entanglement entropy of 
quantum scalar field propagating in flat space-time. In the Appendix, we have shown explicitly 
that at a fixed Lemaitre time, the Hamiltonian of a scalar field in Schwarzschild space-time 
reduces to that in flat space-time. All the results we derive in case of flat space-time can then 
be extended to the black-hole space-time at a fixed Lemaitre time. 

3. Entanglement entropy for scalar fields — Setup 

The Hamiltonian of a free massless scalar field propagating in the Minkowski space-time is 
given by 

H=l I (fx [^2(f) + |Vv9(r)|2] . (14) 



For simplicity, let us assume that the scalar field in confined in a spherical box {12}. The 
cartoon version of the setup is provided in Fig (2). Partial- wave decomposition of the scalar 
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Figure 2. Sketch of the discretization of the scalar field propagating in the flat space-time. 
The shaded region is one we trace over. 

field and its canonical conjugate 



^im{r) = r \ d£l Zim{9,(p)v{r); TTim{r) = r I dTl ZimiO,(j))Tr{f) 



m + il 



leads to the following reduced Hamiltonian: 
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(16) 



where Zim are the real Spherical harmonics. (For details, see Refs. [5, 6, 7].) 

The computation of the entanglement entropy involves three steps: (i) Discretize the 

Hamiltonian, (ii) Choose a quantum state, and (iii) Trace over region 2 (or 1) to obtain the 

density matrix. 

Note that, even if we obtain closed-form expression of the density matrix, it is not possible 

to analytically evaluate the entanglement entropy. Hence, we need to resort to numerical 

methods. In the rest of the section, we discuss these steps in detail. 



3.1. Discretized Hamiltonian 

Discretizing the Hamiltonian (16) in a spherical lattice as described in Fig. (2), we get 
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(17) 
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i,j=i 



where a is the lattice spacing, L (>> a) is the length of the box, N is the number of lattice 
points and the interaction matrix Kij is given by 
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nearest neighbor interaction 
Note that the Hamiltonian (17) is a A^-coupled Harmonic oscillators which satisfy the 



following commutation relations: 
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(18) 



3.2. Choice of quantum state 

The general Eigen state of the A^-coupled harmonic oscillators is 



^(Xl,.. . ,Xn,Xn+l,. ■ .,Xn) 
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where Kd is the diagonal matrix of K obtained by the unitary transformation K = U^KdU, 
U is a unitary matrix and il = U Kjj U . 

It is not possible to obtain a closed form expression for the density matrix for such a general 
state. To make the calculations tractable, we make two choices for the A'"-particle state: 

(i) Vacuum state: In this case, the A^-particle wave function is given by [4] (see also [9, 10]) 
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^q{xi,. .. ,Xrt,X„+l,... ,XAr) 

(ii) 1-particle state: In this case, the A^-particle wave function is given by [5] 

1 1 
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V2 ( a^K^x ) Vo {xi, ..., xn) 



(21) 



where a's are expansion coefficients wliicli satisfy tlie condition a a = 1. For detailed 
discussion see [5, 6]. 
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Figure 3. Tlie entanglement entropy for ground and 1-Particle states for A'^ = 300, n 
100-200,0 = 10-50 



3.3. Density matrix 

For a general Eigen state, the physical operation of trace over region 2 (n < N oscillators) [cf. 

Fig. (2)] is given by 



/n 
Y[ dXi ^(xi . . . XAr)^*(xi ...Xn] X^+1 • • • x'^) 



Pn [Xn+i ■ ■ ■ xn; a;„+i ■ --XN 



For the two special choices which we discussed in the previous subsection, the above 
expression reduces to 

(i) Vacuum state: 



Pas = n /'aC^i'^O 



(22) 



which corresponds to product of n 2-coupled harmonic osciUators with 1 harmonic 
oscihator traced as discussed in Sec. (2). 

1 -Particle state: 
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where 
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(24) 



As mentioned earher, it is not possible to evaluate the entanglement entropy analytically 
and we need to resort to numerical computations. We use Matlab to evaluate the entropy and 
the numerical error in our computation is less than 0.1%. 

In Fig. (3), we have plotted the entanglement entropy for the ground state and 1-particle 
state for N = 300, n = 100 — 200, o = 10 — 50. We see that for the 1-particle state, the entropy 
scales as A^, with lesser ^u for higher o. It is less than unity for any o > i. e., more the 
excitation, larger is deviation from area law The area law does not seem to hold! This implies 
that the entanglement entropy depends on the choice of the quantum state of the field. 



4. Po^ver law corrections to area 

Given the above results, one may draw two distinct conclusions: first that entanglement entropy 
is not robust and reject it as a possible source of black-hole entropy. Second, since entanglement 
entropy for excited state scales as a lower power of area it is plausible that when a generic state 
(consisting of a superposition of ground state and excited state) is considered, corrections to 
the Bekenstein-Hawking entropy will emerge. In order to determine which one is correct, it 
is imperative to investigate various generalizations of the scenarios considered in Ref. [5] . To 
this end, in this section we calculate the entanglement entropy of the mixed superposition of 
vacuum and 1-particle state 



^ = co^o + ci^i 



where 



Fol +|ci| 



1. 



(25) 



Following the procedure discussed in the previous section, it is possible to obtain the 
entanglement entropy of the superposed state. In Fig. (4) we have plotted the relative entropy 
[i.e ratio of the superposed and ground state entanglement entropy] for different values of cq and 
ci. From Fig. (4) it is easy to see that for large values of the horizon radius the entanglement 
entropy of the superposed state approaches the ground state entanglement entropy. In other 
words, for large horizon area, the entanglement entropy of the superposed state scales as area 
of the horizon. 
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Figure 4. Plots of the relative entropy with the area for A^ = 300, n = 100 — 200 and o — 30, 50 
where (i) MSeq corresponds to |co| = |ci| = 0.7 and (ii) MShi corresponds to |co| = 0.5, 
|ci| =0.87. 
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Figure 5. The best fit the relative entropy vs A for (i) MSeq corresponds to |co| = |ci| = 0.7 
and (ii) MShi corresponds to |co| = 0.5, |ci| = 0.87 



In order to know the behavior of the entanglement entropy for the small horizon limit, in Fig. 
(5), we have obtained the best numerical fit for the relative entropy. From these, we see that 
the entanglement entropy of the superposed state is given by: 



S = c,A[l + ^) 



(26) 



where /? > and ci, C2 are constants. This is the main result of this paper and we would like to 
stress the following points: (i) The second term in the above expression may be regarded as a 



power law correction to the area law, resulting from entanglement, when the wave- function of 
the field is chosen to be a superposition of ground and excited states, (ii) It is important to note 
that the correction term falls off rapidly with A (due to the negative exponent) and in the large 
area limit (A ^ 1) the area-law is recovered. This lends further credence to entanglement as 
a possible source of black hole entropy, (iii) The correction term is more significant for greater 
excited state-ground state mixing proportion ci. For detailed discussion see Refs. [7, 8]. 

5. Conclusions 

In this work, we have obtained power-law corrections to the Bekenstein-Hawking entropy, 
treating the entanglement between scalar field degrees of freedom inside and outside the horizon 
as a viable source of black-hole entropy. We have shown that for small black hole areas the area 
law is violated when the oscillator modes are in a linear superposition of ground and excited 
states. We found that the corrections to the area-law become increasingly significant as the 
proportion of excited states in the superposed state increases. Conversely, for large horizon 
areas, these corrections are relatively small and the area-law is recovered. 

Power law corrections to the Bekenstein- Hawking entropy has been encountered earlier in 
other approaches to black hole entropy. For instance, the Noether charge approach predicts 
a generic power law correction to the Bekenstein-Hawking entropy [11]. For instance, using 
Noether charge as entropy, it was shown for 5-dimensional Gauss-Bonnet gravity 

K + acB 



q5D 



IGttG 



d xyj—g 



R' 



A n n^^b _L f? nabcd 



acB IS a constant 



(27) 



(28) 



that the entropy corresponding to a Killing horizon is [12] 

Although both these approaches lead to power-law corrections, there are couple of crucial 
differences between the two: (i) In the case of Noether charge, the power-law corrections to the 
Bekenstein-Hawking entropy arises due to the higher-derivative terms in the gravity action. 
However, the corrections we have derived is also valid for Einstein-Hilbert action, (ii) The 
Noether charge entropy does not identify which degrees of freedom contribute most to the 
entropy. In our case, we have shown that the maximum entropy contribution for the area-law 
come close to the horizon while significant contribution to the power-law corrections arise from 
the degrees of freedom far from the horizon [7, 8] . 
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Appendix A. Scalar fields in Schwarzschild space-time 

In this section, we generalize the framework of our calculations, so that they are applicable to 
black hole space-times. We start with the Schwarzschild metric with a horizon at r = ro: 
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and define the following coordinate transformation from (t, r) to (r, R) coordinates: 
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such that (A.l) transforms to the following metric, in Lemaitre coordinates: 
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(A.3) 



The Hamiltonian for a free scalar field in the above space-time can be written as ^: 
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Next, choose a fixed Lemaitre time, say r = and perform the following field redefinitions: 
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(A.5) 



Then, at that fixed time, the Hamiltonian (A.4) transforms to: 



1 f°^ 
H{0) = - I dr 
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(A.6) 



Note that this is identical to the Hamiltonian in flat space-time, Eq.(14). Hence, all the analysis 
in Sees. (2,3) will go through for a black-hole space-time in a fixed Lemaitre time. 
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